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(57) A method for the determination of horizontal 
resistivity, vertical resistivity, dip and strike angles of 
anisotropic earth formations surrounding a wellbore in 
which electromagnetic couplings among plural triad 
transmitters and triad receivers are measured. Each triad 
transmitter/receiver consists of coil windings in three 
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Method f r Iterative Determination of Conductivity in Anisotr pic Dipping Formations 

BACKGROUND 

Field of the Invention 

The present invention generally relates to the measurement of electrical characteristics of 
formations surrounding a wellbore. More particularly, the present invention relates to a method 
for detennining horizontal and vertical resistivities in anisotropic formations while accounting 
for the dip and stike angle of the formation. 



Description of the Related Art 

The basic principles and techniques for electromagnetic logging for earth formations are 
well known. Induction logging to determine the resistivity (or its inverse, conductivity) of earth 
formations adjacent a borehole, for example, has long been a standard and important technique 
in the search for and recovery of subterranean petroleum deposits. In brief, the measurements 
are made by inducing eddy currents to flow in the formations in response to an AC transmitter 
signal, and then measuring the appropriate characteristics of a receiver signal generated by the 
formation eddy currents. The formation properties identified by these signals are then recorded 
in a log at the surface as a function of the depth of the tool in the borehole. 

It is well known that subterranean formations surrounding an earth borehole may be 
anisotropic with regard to the conduction of electrical currents. The phenomenon of electrical 
anisotropy is generally a consequence of either microscopic or macroscopic geometry, or a 
combination thereof, as follows. 



in many sedimentary strati., electrical eurrent flows more easily in a direction parallel to 
,h« bedding planes, as opposed ,0 a direction Perpendicular to the bedding planes. One reason is 
.hat a great number of mineral crystals possess a flat or elongated shape (e.g., mica or kaolin). 
A, the time they were laid down, mey naturally took on an orientation paralle. to .he plane of 
sedimentation. The interstices in the formations are, therefore, generally parallel to tire bedding 
plane, and the current is able to easily travel along these interstices which often contain 
electrically conductive mineralized water. Such electrical anisotropy, sometimes calfed 
microscopic anisotropy, is observed mostly in shales. 

Subterranean formation are often made up of a series of relatively thin beds having 
different lithological characteristics and, therefore different resistivities. In well logging systems, 
the distances between the electrodes or atuennas are great enough that the volume involved in a 
.neasuremen. may include several such thin beds. When individual layers are neither de.inen.ed 
nor resolved by a togging too!, tire fool responds ,0 the formation as if i, were a macroscopicafly 
anisotropic formation. A thinly laminated sand/shale sequence is a particularly important 
example of a macroscopically anisotropic formation. 

If a sample is cut from a subterranean formation, the resistivity of tire sample measured 
with current flowing parafle. to the bedding planes is called the transverse or horizon*! 
resistivity p.. The inverse of p» is the horizontal conductivity cr„. The resistivity of the sample 
n^asured with a current flowing perpendicular to the bedding plane is called the longitudinal or 
vertienl resistivity, Pv, and its inverse the vertical conductivity <r v . The anisotropy coefficient X is 

defined as: X = ^ja h /o v . 

In situations where the borehole intersects the formation subaantially perpendicular ,0 
the bedding planes, conventional induction and propagation well logging tools are sensitive 



almost exclusively to the horizontal component of the formation resistivity. When the borehole 
intersects the bedding planes at an angle (a deviated borehole) the tool readings contain an 
influence from the vertical and horizontal resistivities. This is particularly true when the angle 
between the borehole and the normal to the bedding places is large, such as in directional or 
horizontal drilling, where angles near 90° are commonly encountered. In these situations, the 
influence of vertical resistivity can cause discrepancies between measurements taken in the same 
formation in nearby vertical wells, thereby preventing a useful comparison of these 
measurements. In addition, since reservoir evaluation is typically based on data obtained from 
vertical wells, the use of data from wells drilled at high angles may produce erroneous estimates 
of formation reserve, producibffity, etc. if proper account is not taken of the anisotropy effect. 

There have been proposed a number of methods to determine vertical and horizontal 
resistivity near a deviated borehole. Hagiwara (US Patent 5,966,013) disclosed a method of 
determining certain anisotropic properties of formation using propagation tool without a priori 
knowledge of the dip angle. In US patent 5,886,526, Wu described a method of detennining 
anisotropic properties of anisotropic earth formations using multi-spacing induction tool with 
assumed functional dependence between dielectric constants of the formation and its horizontal 
and vertical resistivity. Gupta et al. (US 5, 999, 883) utilized a triad induction tool to arrive at an 
approximate initial guesses for the anisotropic formation parameters. Moran and Gianzero 
(Geophysics, Vol. 44, P. 1266, 1979) proposed using a tri-axial tool of zero spacing to 
determine dip angle. Later the spacing was extended to finite size by Gianzero et al. (US Patent 
5, 1 15, 198) using a pulsed induction tool. The above references are hereby incorporated herein 
by reference. 



These attempts to determine vertical and horizontal resistivity around a deviated 
borehole have thus far not provided sufficient accuracy for formations having a high degree of 
anisotropy. A new technique is therefore needed. 

SUMMARY OF THE INVENTION 
The above-described problems are in large part addressed by an iterative method for 
determining electrical conductivity in an anisotropic dipping formation. The iterative method 
corrects for the skin effect to high orders while detennining all relevant formation parameters. 
This method may be applied to a tri-axial induction sonde operating in continuous wave (CW) 
mode. In one embodiment, the method includes (1) measuring a magnetic coupling between 
transmitter coils and receiver coils of a tool in a borehole traversing the formation; (2) obtaining 
from the measured coupling a strike angle between the tool and the formation; (3) obtaining 
from the measured coupling an initial dip angle between the tool and the formation; (4) obtaining 
from the measured coupling an initial anisotropic factor of the formation; (5) obtaining from the 
measured coupling an initial horizontal conductivity of the formation; (6) detennining an 
iterative anisotropic factor from the measured coupling, the strike angle, the latest dip angle, and 
the latest anistropic factor; (7) determining an iterative horizontal conductivity from the 
measured coupling, the strike angle, the latest iterative anisotropic factor, and the latest dip 
angle; and (8) determining an iterative dip angle from the measured coupling, the latest iterative 
anisotropic fector, and the latest iterative horizontal conductivity. The steps of detennining an 
iterative anisotropic factor, determing an iterative horizontal conductivity, and determining an 
iterative dip angle are preferably repeated a number of times that minimizes an overall residual 
error. 



The disclosed method may provide the following advantages in determining the 
formation parameters of anisotropic earth formations: (1) a priori knowledge of the dip angle is 
unnecessary and can be one of the outputs of the method; (2) no assumed relationship between 
formation resistivity and dielectric constant is necessary; (3) complex electronics for pulsing the 
transmitter coils may be eliminated since this method is applicable to a triad induction sonde 
running in CW mode; (4) preliminary results indicate that the disclosed method yields more 
accurate estimates of all electrically relevant formation parameters in the earth formation. 

BRIEF DESCRIPTION OF THE DRAWINGS 

A better understanding of the present invention can be obtained when the following 
detailed description of the preferred embodiment is considered in conjunction with the following 
drawings, in which: 

Fig. 1 shows the coil configuration of a triazial induction tool; 

Fig. 2 demonstrates a rotational transformation definition; 

Fig. 3 shows a flow diagram for the disclosed method of determining formation 

parameters in a dipping anisotropic earth formation; 

Fig. 4 shows a graph used to determine the optimal number of iterations that minimize a 
dip-angle error function; and 

Fig. 5 compares the results of the disclosed method to results obtained with an existing 
method. 

While the invention is susceptible to various modifications and alternative forms, specific 
embodiments thereof are shown by way of example in the drawings and will herein be described 
in detail. It should be understood, however, that the drawings and detailed description thereto 



are not intended to limit the invention to the particular form disclosed, but on the contrary, the 
intention is to cover all modifications, equivalents and alternatives falling within the spirit and 
scope of the present invention as defined by the appended claims. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Terminology 

It is noted that the terms horizontal and vertical as used herein are defined to be those 
directions parallel to and perpendicular to the bedding plane, respectively. 
Tool Configuration 

Turning now to the figures, Fig. 1 shows a conceptual sketch of a coil arrangement for a 
downhole induction tool. A triad of transmitter coils T x , T y and T z , each oriented along a 
respective axis, is provided, as is a similarly oriented triad of receiver coils R x , R y and R z ,. For 
clarity, it is assumed that the three coils in each triad represent actual coils oriented in mutually 
perpendicular directions, with the z-axis corresponding to the long axis of the tool. However, it 
is noted that this coil arrangement can be "synthesized" by performing a suitable transformation 
on differently oriented triads. Such transformations are described in depth in U.S. Patent 
Application Serial No. 09/255,621 entitled "Directional Resistivity Measurements for Azimutal 
Proximity Detection of Bed Boundaries" and filed February 22, 1999 by T. Hagiwara and H. 
Song, which is hereby incorporated herein by reference. 

Each of the coils in the transmitter triad is parallel to the corresponding coil in the 
receiver triad and is spaced away from the corresponding coil in the z-axis direction. The 
distance between the corresponding coils is labeled L. It is noted that the downhole tool may 
have additional transmitter or receiver triads to provide multiple transmitter-receiver spacing L 



values. Such configurations may advantageously provide increased accuracy or additional detail 
useful for analyzing the formation structure. 
System Model 

Generally, a formation model is used to interpret the tool measurements. The model used 
herein is a unixial anisotropy model. This model assumes that the formation is isotropic in the 
horizontal direction (parallel to the bedding plane) and anisotropic in the vertical direction 
(perpendicular to the bedding plane). Setting up a formation coordinate system having the z-axis 
perpendicular to the bedding plane and the x- and y-axes parallel to the bedding plane allows a 
conductivity tensor to be expressed as: 



cr h 0 0 
0 <r A 0 
0 0 <t 



(1) 



The axes of the formation coordinate system typically do not correspond to the axes of 
the tool coordinate system. However, a rotational transformation from one to the other can be 
defined. Fig. 2 shows a transformation from the tool coordinate system to the formation 
coordinate system. The tool coordinate system (x,y,z) is first rotated about the z-axis by an 
angle (5, hereafter termed the strike angle. The intermediate coordinate system (x\y',z'=z) thus 
formed is then rotated about the y* axis by an angle a, hereafter termed the dip angle to obtain 
the formation coordinate system (x",y"=y\z")- 

Any vector v" in the formation coordinate system can be expressed in the tool 
coordinate system as: 

v = Rv" (2) 
where the rotational transform matrix is: 



R = 



cos a cos P cos a sin P —sin a 

-sin a cosp 0 
sin a cos p sin a sin p cos a 



(3) 



Derivation of Equations 

Now that the rotational transformation has been defined, attention is directed to the 
induction tool measurements. When a voltage is applied to one of the transmitter coils, a 
changing magnetic field is produced. The magnetic field interacts with the formation to induce a 
voltage in the receiver coils. Each of the three transmitter coils is excited in turn, and the voltage 
produced at each of the three receiver coils is measured. The nine measured voltages indicate 
the magnetic coupling between the transmitter-receiver triad pair. The equations for the 
measured signals will be derived and manipulated to solve for the strike angle P, the dip angle a, 
the horizontal conductivity Oh, and the vertical anisotropy X. 

In the most general case according to Moran and Gianzero (Geophysics, Vol. 44, P. 
1266, 1979), the magnetic field H in the receiver coils can be represented as a coupling matrix C 
in the form: 



where H x , Hy, H- and M x , M y , A4, are the field components at the receivers and magnetic 
moment components at the transmitters, respectively. (Magnetic moment is calculated 
Mt=AtNtIt, where A T is the transmitter area, N T is the number of turns in the transmitter coil, 
and I T is the transmitter current. The direction of the magnetic moment is perpendicular to the 
plane of the coil). If the coupling matrix is specified in terms of the formation coordinate system, 




C C C 

xx xy xz 

c c c 

^yx ^yy ^zz 

c B C B 



M x 
M y 



(4) 



the measured magnetic field strengths in the receiver coils are obtained by (Moran and Gianzero, 
Geophysics, Vol. 44, P. 1266, 1979): 

h = (r-'cr)m=cm ( 5 > 

where H and M are measured in the sonde coordinate system 

Assuming that the tool is oriented so that the strike angle p is 0, it can be shown that for 
the uniaxial anisotropy model the full coupling matrix C» at the receiver coils (x =0 =y, z =L) 
simplifies to 



C'= 



c„ o c a 

0 C„ 0 

CL 0 C„ 



(6) 



The theoretical values of the coupling matrix elements are (Q- Qi): 

1 kit cos 2 a-e^ ~tt,4 ^T^ Vg« 2 «- 1 )0-- <fc ^H (?) 

^ sin a y J 



C = 



Anlf 



C *y= C >*= 0 



_ -sin a cos a 
C = s e 



AnU 



^\-ik h L)-klL 2 ] 



C » ~ Aid? 



an 

XA 1 



e iKL +ih h L\ 



sin 2 a 



(l-i* A I)e^j 



C = — Ure^ fex 2 sin 2 a + (1 - /A: A I)(3cos 2 a - 1)] 



(8) 
(9) 

(10) 
(11) 



where 



k h = yji(ay.a h = horizontal wave number 

© = 2nf = angular frequency 

p. = p. 0 = 4% x 1 0" 7 henry I m = magnetic permeability 
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X - ^ja h /c v = anisotropy coefficient 

A = Vsin 2 a + A, 2 cos 2 a/x, = anisotropic factor 
In terms of the elements of the coupling matrix C\ the six independent measurements for 
all the possible couplings between all transmitter-receiver pairs are expressed as (T,Rj=TjRj): 



^r=y{[(C«+C„) + (C„ - C„ ) cos2a + 2C„ sin 2a]cos 2 p + 2C^ sin 2 p } 



(12) 



T y R y = Y { [ ( c « + C s ) + (C„ - C„) cos2a + 2C„ sin 2a]sin 2 P + 2C„, cos 2 p } (13) 



M r 

r A = -=ri (C„ + C s ) + (C a - C„)cos 2a - 2C„ sin 2a] 



(14) 



r A=x( (C « +C = ) + (C «- C -) cos2a + 2C « sin2a - 2 ^]sin2p (15) 



A/r 

^ =—[(C 22 -C^)sin2a+2C„ cos2a]cosp 



(16) 



Mr 

T z R y =—[(C zz -C„)sin2a+2C l2 cos2a]sinp (17) 

These measurements are made by taking the ratio of the transmit and receive voltage signals, 
e.g. T x Ry = K Vr/Vtx, where K is a real-valued calibration constant theoretically equal to 
A r N r I T A n N R (a)/j) 2 /4xL , where A R is the area of the receive coil, and N R is the number of turns 
of the receive coil. 

Explicitly solving the last four of the above equations results in the following expressions 
for the measured cross-coupling fields: 



TR - Msm 2 P 



1 y 



( J^IA \ 



X 2 A 



-ik h L(e^-e^ 1 + ct f a ) 



sin a 



(18-a) 
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T 2 R X = Mc0 ? sin 2 2a ik h L(e*>» -e*+) 
y 8nZ, 3 sin 2 a * v ' 



(18-b) 
(18-c) 
(18-d) 



To make practical use of the above equations, the real component is ignored and the 
imaginary (reactive) component is simplified by finding the limit as the transmitter-receiver 
spacing approaches zero, i.e., L -> 0. Doing this simplifies the reactive components of the 
measured signal equations (18-a,b,c) to: 



( ) ^ Msm2 P sina , 



(rA) _^M^s|2a (1 _ x2) 



( ) ^ MsinPsin2a 2 
v 2 y,x %itiZh\ 



(19-a) 
(19-b) 
(19-c) 



where 5 A = ^2/(a\xcs h is the skin depth associated with horizontal conductivity. From these 
equations, one arrives at the practical equations for the determination of dip and strike angles: 



/? = tan 



a ol = tan" 



1 imii 


2 




2 











(20) 



(21) 
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It is noted that the strike angle P thus obtained is exact while the dip angle a is only an 
approximation because Equations (19a-c) are valid only in the zero-spacing limit. The subscript 
a\ denotes that this is the first approximation of the apparent dip angle. 

With the strike angle P and the estimated dip angle a, estimates of the horizontal 
conductivity an and anisotropy factor A can be obtained via the following observation. 

When a power series expansion is used for the exponential terms in equation (18-b), the 
first terms yield the following expressions for the real part (T z R x )r and imaginary part (T Z R X ) X : 

M cos p sin 



8tiZ 3 sin 



M cos P sin 2a 



8jcL sin a 



-2{A-\) 



+ (A 2 -l) 



f4 



(22) 



(23) 



Taking advantage of the feet that the second term in (T Z R X ) X is identical to (T z R x )r, an equation 
that is skin-effect corrected to the first order can be written as: 



<W. -<W, . " ""Pf 3 " (A -1) 

4kL sin a 



Similarly with equation (18-d), for T Z R Z one has: 



'4 



(24) 



M 
4nL 3 



(25) 



Equations (24) and (25) can be solved to obtain A and ov 
1 



A*- 



1 + 



(T Z R X ) X -(T,R X ) R tana 



a, « 



(T 2 R 2 ) x -(T t R z ) R cosp 



(26) 



(27) 



- 12- 



Substituting the strike angle p and the first estimate of the dip angle a„, from Equations (20) and 
(21) yields the first estimates of the anisotropic factor A a i and the horizontal conductivity o AI , 
where the subscript a\ denotes the quantities are the first estimates of the apparent values: 

1 



n (T Z RJ X -(T Z R X ) R tang,, 
i T ,K)x-(7' 2 R I ) R * cos/? 

M<o\x A a , 



(28) 



(29) 



It is noted that the leading fraction may be replaced by a calibration coefficient K, for the tool. 

Now that initial estimates have been obtained, the estimates can be iteratively improved. 
Examining the series expansion of all the measured fields reveals that 



_ M cos/?s in(2a) 



sin 2 a 



(^-1) \±- 



rr p \ - M <*>sy?sin(2a) 



sura 



(30) 



(31) 



Because the higher order terms do not cancel each other, Equations (24) and (25) are no longer 
correct. To remedy this problem, the extra skin-effect contribution may be subtracted from the 
left-hand side. Namely, the left-hand side of Equation (24) may be replaced by {T Z RJ R -(T Z R X ) X - 
T a where T a is the higher order correction obtained from (30): 



r M cosPsin2q al 



8tiL sin 2 a 



at 



(32) 
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Similarly, the correction for (TzRJr -(T t R z )x can be derived: 



r = 



M 
4nL 3 



-L(l + 2A al 3 ) 



f 



0 + 5Aj)( L 



15 



(33) 



It is noted that the leading fraction in (32) and (33) may be replaced with a calibration 
coefficient K 2 and 2K2, respectively, for the tool. With these corrections, one gets better 
approximations for A and at,'. 

1 



A ° 2 (TM x -(T z R x ) R -r a tana c 



1 + 



°H7 = 



(J,R t ) x -(T z R x ) R -T a cosp 
_ 4*L [(T Z R X ) X -(T Z R Z ) R ] 



Mcofj 



(34) 



(35) 



1 a2 



The same "extend to a higher order" procedure can be applied to the other components 
of the measured fields in Equation (25). The end result is a more accurate equation for the dip 
angle: 



tana,, = . 



(t z rX 



(T x R y ) R ~(K">* 



(36) 



It is noted that here the real part of the magnetic field (which is identical to the imaginary part 
of the measured voltages at the coils other than a constant factor) is used, hence the subscript R. 
The correction terms can be directly obtained from the 6 th order expansions of the corresponding 
coupling fields. The correction terms are: 



/a x Af cospsin(2a g ) 



SnL sin cl 



M sinpsin(2a fl ) 



8tiL 3 



sin a. 



(37-a) 



(37-b) 
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Msm2g 







4 






f L ] 















(37-c) 



where 



= Apparent dip angle 

Apparent horizontal skin-depth 



Oha = Apparent horizontal conductivity 

p(pLa) „ (4,* -1) 0+cos 2 qJ t A a (A a 2 -cos 2 cc J 1 
6 sin 2 a a 2 sin 2 a„ 2 



g(aJ _ K 4 -D(l + cos 2 aJ | V (4, 2 -cos 2 aJ 1 
24 sin 2 a /7 6 sin 2 a„ 6 



Property Determination Method 

The above derivation has provided: 

Eqn. (20) for the determination of the strike angle 0, 

Eqn. (21) for a first estimate of the dip angle a, 

Eqn. (28) for a first estimate of the anisotropic factor A 9 

Eqn. (29) for a first estimate of the horizontal conductivity a, 

Eqn. (34) for an iterative estimate of the anisotropic factor A, 

Eqn. (35) for an iterative estimate of the horizontal conductivity a, and 

Eqa (36) for an iterative estimate of the dip angle a. 

Eqn. (36) includes correction terms that are specified by Eqn. (37). The property determination 

approach built into Eqns (34)-(37) is believed to be unkown to the art. 

Fig. 3 shows an iterative method that uses the above equations to determine the dip and 

strike angles, the horizontal conductivity, and the anisotropy coefficient (calculated from the 

anisotropic fector). In block 302, the transmitter coils are activated and the voltages produced in 

the receiver coils are measured. The measured voltages are processed to determine the magnetic 
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couplings between the coils. In block 304, the strike angle p is calculated per Eqn (20) and the 
first estimate of the dip angle a A i is calculated per Eqn (21). In block 306, the first estimate of 
the anisotropic factor A A is calculated per Eqn (28) and the first estimate of the horizontal 
conductivity is calculated per Eqn (29). 

Blocks 308 through 3 14 form a loop having a loop index J. Loop index J is initialized to 
1 for the first iteration of the loop, and incremented by 1 for subsequent iterations of the loop. 

In block 308, an iterative estimate of the anisotropic factor Aj is calculated per Eqn (34) 
and an iterative estimate of the horizontal conductivity ow is calculated per Eqn (35). In block 
310, coupling correction terms are caculated per Eqns (37-a,b,c). In block 312, an iterative 
estimate of the dip angle a^n is calculated per Eqn (36) using the coupling correction terms 
found in block 310. In block 314, a test is made to determine if the optimum number of 
iterations has been performed, and if more iterations are needed, control returns to block 308. 
Otherwise, in block 316, the anisotropy coefficient X is calculated from the latest estimates of 
the dip angle and anisotropic feictor (refer to the definition of the anisotripic factor in Eqn (1 1)). 
The anisotropy coefficient X is then output, along with the strike angle 0, and the best estimates 
of the dip angle a and horizontal conductivity o h . 
Analysis of Method 

The effectiveness of the method is demonstrated in Table I where all apparent formation 
parameters are listed for each successive iteration. The true formation parameters are: dip angle 
a = 45°, strike angle p = 60°, anisotropy coefficient X = 2, and horizontal conductivity o h = 5 
S/m. The sonde is a triad 2C40 (notation means 2 coils spaced 40 inches apart) operating at 8 
kHz. 

Table I: Effect of iteration number for triad 2C40 in a dipping anisotropic formation 



- 16- 



ILCrdllOIl 1MO. 


/\pparent uip cto 


Apparent A, 


Apparent 


A 

u 


34. /2 


A *J1 

4.71 


9.234 


1 
1 


3y.21 


1.916 


4. /JZ 


2 


yf 1 "7/1 ° 


1 A/'O 

1.968 


4.590 


3 


43.22 


1 .992 


4.975 


A 


44.1 1 




3.U29 


5 


44.65° 


2.013 


5 060 


6 


44.99° 


2.017 


5.083 


7 


45.19° 


2.020 


5.095 


8 


45.32° 


2.021 


5.104 


9 


45.40° 


2.022 


5.109 


True Values 


45.00° 
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A close examination of Table I reveals that there exists an optimal iteration number for 
each formation parameter beyond which the accuracy of apparent formation parameter actually 
decreases. For a given triad sonde, the optimal number of iterations for a given parameter can be 
calculated with the following method. For the three formation parameters, their corresponding 
error functions are defined as: 



y 
) 



(38-a) 



(38-b) 



e *» =ZZZ 

a a h X 



(38-c) 



The summations are performed over the range of expected formation dip angles, conductivities, 
and anisotropy coefficients. This error function is consequently indicative of the overall residual 
error. 

Given a configuration of a triad sonde, the functional dependence of the error function 
on the number of iterations can be examined. Fig. 4 shows the apparent dip angle examination 
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results for four triad-pair configurations having different transmitter-receiver spacings. The error 
function is calculated over the following ranges of formation parameters: 
a = (5, 15, 30, 45, 60, 75, 85°) 

a„ = (0.001, 0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1, 2, 5, 10 S/m) 
X=(l.l, 1.414, 2,4,6,8, 10). 
The optimal iteration number is defined as the iteration number at which the error is rninimal. 
This can similarly be done for the anisotropy coefficient X and horizontal conductivity 07, to 
obtain their respective optimal iteration numbers. 

With the optimal iteration number determined, the advantage of the current invention in 
the dip angle determination over an existing method (Moran and Gianzero, Geophysics, Vol. 
44, P. 1266, 1979) can be demonstrated. Fig. 5 plots the calculated dip angle measured in 
dipping formations having different horizontal conductivities. The true dip angle of the 
formation is 75 degrees. It is noted that the original method proposed by Moran and Gianzero 
has large errors in conductive formations due to the skin effect, while the disclosed method 
yields a calculated dip angle close to the true dip of 75 degrees. It is further emphasized that in 
addition to a more accurate dip angle, the current method also yields the other formation 
parameters illustrated in Table I. 

Numerous variations and modifications will become apparent to those skilled in the art 
once the above disclosure is folly appreciated. It is intended that the following claims be 
interpreted to embrace all such variations and modifications. 
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CLAIMS 

WHAT IS CLAIMED IS: 

1. A method for detennining conductivity in a formation, wherein the method comprises: 

measuring a magnetic coupling between transmitter coils and receiver coils of a tool in a 

borehole traversing the formation; 
obtaining from the measured coupling a strike angle between the tool and the formation; 
obtaining from the measured coupling an initial dip angle between the tool and the 

formation; 

obtaining from the measured coupling an initial anisotropic factor of the formation; 
obtaining from the measured coupling an initial horizontal conductivity of the formation; 
determining an iterative anisotropic factor from the measured coupling, the strike angle, 

the latest dip angle, and the latest anistropic factor; 
determining an iterative horizontal conductivity from the measured coupling, the strike 

angle, the latest iterative anisotropic factor, and the latest dip angle; and 
determining an iterative dip angle from the measured coupling, the latest iterative 

anisotropic factor, and the latest iterative horizontal conductivity. 

2. The method of claim 1, further comprising: 

repeating the steps of determining an iterative anisotropic factor, determing an iterative 
horizontal conductivity, and determining an iterative dip angle. 

3. The method of claim 2, wherein said repeating is performed a number of times that minimizes 
an overall residual error. 
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4. The method of claim 1 , wherein the strike angle p corresponds to 

wherein (T-R y ) x is the reactive component of the coupling T,R y between a transmitter T z oriented 
along a z-axis and a receiver Ry oriented along a y-axis, and (T Z R X ) X is the reactive component of 
the coupling T Z R X between transmitter T 2 and a receiver R x oriented along an x-axis. 

5. The method of claim 1, wherein the initial dip angle ai corresponds to 



frfc^ii 


2 


"fc*J,T 




+ 





wherein (T x R y ) x is the reactive component of the coupling T x R y between a transmitter T x 
oriented along an x-axis and a receiver R y oriented along a y-axis, (T Z R X ) X is the reactive 
component of the coupling T Z R X between a transmitter T z oriented along a z-axis and a receiver 
R x oriented along the x-axis, and (T z R y ) x is the reactive component of the coupling T z R y between 
transmitter T z and receiver Ry. 

6. The method of claim 1, wherein the initial anisotropic factor A, corresponds to 

A = ! 

1 U (T Z RJ X -(T:R X ) R tana/ 

(T Z R 2 ) X -{T 2 R 2 ) R cosp 
wherein (T Z R X ) X and (T Z R X ) R are the imaginary and real components, respectively, of the coupling 
T Z R X between a transmitter T z oriented along a z-axis and a receiver R x oriented along an x-axis, 
(T Z R Z ) X and (T Z R Z ) R are the imaginary and real components, respectively, of the coupling T Z R Z 
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between transmitter T z and a receiver ^ oriented along the z-axis, ai is the initial dip angle, and 
p is the strike angle. 

7. The method of claim 1, wherein the initial horizontal conductivity o h i corresponds to 

°h\ — ^1 - > 

A 

wherein (T Z R X ) X and (T Z R X ) K are the imaginary and real components, respectively, of the coupling 
T Z R X between a transmitter T z oriented along a z-axis and a receiver Rx oriented along an x-axis, 
Ai is the initial anisotropic factor, and Kj is a predetermined function of transmitter signal 
voltage and frequency. 

8. The method of claim 1, wherein the iterative anisotropic factor Aj+i corresponds to 

1 



4+i — 



u (T 2 R JC ) x -(T z R x ) R -r xx tana, 9 
(T z R 2 ) x -(T z R 2 ) R -r B cos/3 



wherein (r*R x )x and (T Z R X ) R are the imaginary and real components, respectively, of the coupling 
T Z R X between a transmitter T z oriented along a z-axis and a receiver R x oriented along an x-axis, 
(T Z R 2 ) X and (T Z R Z ) K are the imaginary and real components, respectively, of the coupling T Z R Z 
between transmitter T z and a receiver R 2 oriented along the z-axis, a { is the latest dip angle, ji is 
the strike angle, T^ is a first skin-effect correction, and T^ is a second skin effect correction. 

9. The method of claim 8, wherein the first skin-effect correction r„ corresponds to 



„ cos j£ sin 2a; 

1 zx - A 2 ~j 



sin a i 



fd-4 3 ) 





5 


(*) 









and the second skin-effect correction corresponds to 
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r =2#, 



wherein L is a distance between the transmitter and receiver, Ai is the latest anistropic factor, 



& hi — ^f^^hi *s toe latest skin depth, a hi is the latest horizontal conductivity, and K 2 is a 
predetermined function of transmitter signal voltage. 

10. The method of claim 1, wherein the iterative horizontal conductivity Oh(i+i) corresponds to 

_ J , | to-(wJ ) 



'*(/+!) 



A. 



i+1 



wherein (T : R X ) X and (T Z R X ) R are the imaginary and real components, respectively, of the coupling 
T Z R X between a transmitter T 2 oriented along a z-axis and a receiver R x oriented along an x-axis, 
Ahi is the latest anisotropic factor, L is a distance between the transmitter and receiver, and Kj 
is a predetermined function of transmitter signal voltage and frequency. 

1 1 . The method of claim 1, wherein the iterative dip angle a i+ i corresponds to 



<**i=tan 



-i 



wherein (T x R y ) R is the real component of the coupling T x R y between a transmitter T x oriented 
along an x-axis and a receiver R y oriented along a y-axis, (T Z R X ) R is the real component of the 
coupling T Z R X between a transmitter T z oriented along a z-axis and a receiver R x oriented along 
the x-axis, (T z R y ) R is the real component of the coupling T z R y between transmitter T z and 
receiver Ry, (A xy ) R is a first correction term, (A^r is a second correction term, and (A zy ) R is a 
third correction term 
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12. The method of claim 11, wherein the first correction term (AxyV corresponds to 

\ 4 



= Msin2/7 



-Pipe,) 



— T 



6 sin a, 2 sin^a, 2 



24 



sin 2 of; 



sin 2 a, 



wherein the second correction term (Azx)r corresponds to 



M cos^sin(2a j ) 



8tcL sin a ; 



v 



+ i(4 +1 4 -l) 



and wherein the third correction term (A zy )R corresponds to 



otzL sm a, 



3 v<WoJ 6 



wherein a t is the latest dip angle, 0 is the strike angle, L is a distance between the transmitter and 



receiver, A i+ i is the latest anistropic factor, £ A(i+1) = ^2/a>//cr M/+l) is the latest skin depth, o h (i+i) is 

the latest horizontal conductivity, and K 2 is a predetermined function of transmitter signal 
voltage. 



13. The method of claim 1, wherein the transmitter coils consist of a triad of mutually 
orthogonal transmitters. 

14. The method of claim 13, wherein the receiver coils consist of a triad of mutually orthogonal 



receivers. 
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15. The method of claim 1, wherein said measuring includes exciting each transmitter coil in turn 
and measuring in-phase and quadrature phase voltage signals induced in each of the receiver 
coils by each of the transmitter coils. 

1 6. A method for detennining conductivity in a formation, wherein the method comprises: 

receiving magnetic coupling measurements from an induction tool; 

obtaining from the magnetic coupling measurements an initial horizontal conductivity of 
the formation; 

determining an iterative horizontal conductivity from the magnetic coupling 
measurements and the initial horizontal conductivity of the formation. 

1 7. The method of claim 16, wherein the initial horizontal conductivity o hl corresponds to 

A 

wherein K, is a predetermined function of transmitter signal voltage and frequency, (T Z R X ) X and 
(7^ are the imaginary and real components, respectively, of the magnetic coupling 
measurement T Z R X between a transmitter T 2 oriented along a z-axis and a receiver R x oriented 
along an x-axis, and A, is the initial anisotropic factor corresponding to 

1 



A = 



u (T,R x ) x -(T t R r ) B tana,' 
(T 2 R,) X -(T 2 R 2 ) R ' cos p 



wherein wherein the strike angle p corresponds to 



= tan 

and the initial dip angle a, corresponds to 
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a, = tan" 



TfcU 1 


2 






+ 





(7W? z )x and (T Z R Z ) R are the imaginary and real components, respectively, of the coupling T Z R Z 
between transmitter T 2 and a receiver R* oriented along the z-axis, wherein (T x R y ) x is the 
reactive component of the coupling T x R y between a transmitter T x oriented along the x-axis and 
a receiver R y oriented along a y-axis, wherein (T z R y ) x is the reactive component of the coupling 
T z R y between transmitter T 2 and receiver R y . 

18. The method of claim 16, wherein the iterative horizontal conductivity o h(1+1 ) corresponds to 

; > 



'*(/+!) 



wherein K, is a predetermined function of transmitter signal voltage and frequency, wherein 
(T Z R X ) X and (T 2 R X ) K are the imaginary and real components, respectively, of the coupling T Z R X 
between a transmitter T, oriented along a z-axis and a receiver R x oriented along an x-axis, and 
A i+ i corresponds to 



i , (r.*,Xr-(rA),-r, tang, ' 

iT z R z ) x -{T z R z ) R -T B cos/? 
wherein (T Z R Z ) X and (T Z R Z ) R are the imaginary and real components, respectively, of the coupling 
T Z R Z between transmitter T 2 and a receiver R 2 oriented along the z-axis, r a is a first skin-effect 
correction, 1^ is a second skin effect correction, the strike angle 0 corresponds to 



p = tan 



and the dip angle a, corresponds to 
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a, = tan 



wherein (T 2 R y )x and (r^) R are the imaginary and real components, respectively, of the coupling 
T z R y between transmitter T 2 and a receiver R y oriented along a y-axis, (T x R y ) R is the real 
component of the coupling T x R y between a transmitter T x oriented along the x-axis and receiver 
Ry, (A X y)R is a first correction term, (A zx )r is a second correction term, and (Azy) R is a third 
correction term. 

19. The method of claim 18, wherein the first skin-effect correction corresponds to 



T- _ ir <x>s/?sin 2a , 
sin a k 



and the second skin-effect correction corresponds to 



jO + 24 3 ) 



f L V (3 + 54 2 > r 



wherein L is a distance between the transmitter and receiver, Ai is the latest anistropic factor, 
S hi = ^2/ copier hi is the latest skin depth, a h j is the latest horizontal conductivity, and K2 is a 
predetermined function of transmitter signal voltage. 



20. The method of claim 18, wherein the first correction term (A xy ) R corresponds to 



, K x Msm2B 







4 






f z ] 












UJ. 



with j>( g .) - <4+i 3 ~ D 0 + cos2 g -> 1 4*1 (4 +l 2 ~ cos 2 Of,) 1 ^ 
6 sin 2 a, 2 sin 2 «. 2' 
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24 sin 2 a, 6 sin 2 a, 6 



wherein the second correction term (A^r corresponds to 



/a ^ - M cosygsinC^.,) 



$7iL sin a f _, 



f(i-4 3 ) 



and wherein the third correction term (A zy ) R corresponds to 



M sinygsin(2g,_,) 
8;rf. sin a w 



(>-4 J ) 



1<V 



4(^-1) 



1<W 



wherein cm is the latest dip angle, p is the strike angle, L is a distance between the transmitter and 
receiver, Aj is the latest anistropic factor, 8 hi = ^2/(ojuar hi is the latest skin depth, ow is the latest 



horizontal conductivity, Ka is a predetermined function of transmitter signal voltage. 
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